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1. Introduction 

This paper presents an implementation of a new motion planning algorithm for general 
robot manipulators moving among three-dimensional polyhedral obstacles. The algo¬ 
rithm has a number of advantages: it is simple to implement, it is fast for manipulators 
with few degrees of freedom, it can deal with manipulators having many degrees of free¬ 
dom (including redundant manipulators), and it can deal with cluttered environments 
and non-convex polyhedral obstacles. An examples of a path obtained from an imple¬ 
mentation of the algorithm is shown in Figure 1. 

The ability to automatically plan collision-free motions for a manipulator given 
geometric models of the manipulator and the task is one of the capabilities required to 
achieve task-level robot programming [15]. Task-level programming is one of the principal 
goals of research in robotics. It is the ability to specify the robot motions required to 
achieve a task in terms of task-level commands, such as “Insert pin-A in hole-B”, rather 
than robot-level commands, such as “Move to 0.1,0.35,1.6”. 

The motion-planning problem, in its simplest form, is to find a path from a specified 
starting robot configuration to a specified goal configuration that avoids collisions with a 
known set of stationary obstacles. Note that this problem is significantly different from, 
and quite a bit harder than, the collision detection problem: detecting whether a known 
robot configuration or a path would cause a collision [l, 4]. Motion planning is also 
different from on-line obstacle avoidance : modifying a known robot path so as to avoid 
unforeseen obstacles [6, 9, 10, ll]. 

Although general-purpose task-level programming is still many years away, some 
of the techniques developed for task-level programming are relevant to existing robot 
applications. There is, for example, increasing emphasis among major robot users on 
developing techniques for off-line programming, by human programmers, using CAD 
models of the manipulator and the task. In many of these applications motion plan¬ 
ning plays a central role. Arc welding is a good example; specifying robot paths for 
welding along complex three-dimensional paths is a time-consuming and tedious process. 
The development of practical motion-planning algorithms could reduce significantly the 
programming time for these applications. 

A great deal of research has been devoted to the motion-planning problem within 
the last five to eight years, e.g., [2, 3, 5, 7, 8, 12, 13, 14, 16, 17, 19, 20]. But, few of these 
methods are simple enough and powerful enough to be practical. Practical algorithms 
are particularly scarce for manipulators made up of revolute joints, the most popular 
type of industrial robot. I know of only two previous motion-planning algorithms that 
are both efficient and reasonably general for revolute manipulators with three or more 
degrees of freedom [2, 7]. Brooks’s algorithm [2] has demonstrated impressive results, 
but is fairly complex. Faverjon’s algorithm [7], on the other hand, is appealingly simple. 
The basic approach of the algorithm described here is closely related to the method 
described by Faverjon. Many of the details of the present algorithm, however, especially 
the treatment of three-dimensional constraints and the free space representation, are new 
and more general. 




Figure 1. (a) Sample configurations along a path for all six links of a Unimation Puma manipulator 
obtained using the algorithm described here. The total planning time on a Symbolics 3600 Lisp 
Machine without floating-point hardware was 40 seconds. The incremental cost for planning 
another path in the same environment is approximately 2 seconds. 
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The approach taken in this algorithm is similar to that of [7, 8, 12, 13] in that it 
involves quantizing joint angles. It differs in this respect from exact algorithms such as 
[17, 19]. On the other hand, the quantization approach lends itself readily to efficient 
computer implementation. 

The purpose of this paper is to show that motion planning for general manipulators 
can be both simple and relatively efficient in most practical cases. There is no reason 
why motion planning should be any less practical than computing renderings of three 
dimensional solids in computer graphics. In both cases, there are many simple numerical 
computations that can benefit from hardware support. In fact, it is worth noting that 
in the examples in Figure 1 it took longer to compute the hidden-surface displays in the 
figures than to compute the paths. 

2. The Basic Approach: Slice Projection 

The configuration of a moving object is any set of parameters that completely specify the 
position of every point on the object. Configuration space (C-space) is the space of con¬ 
figurations of a moving object. The set of joint angles of a robot manipulator constitute 
a configuration. Therefore, a robot’s joint space is a configuration space. The cartesian 
parameters of the robot’s end effector, on the other hand, do not usually constitute a 
configuration because of the multiplicity of solutions to a robot’s inverse kinematics. It 
is possible to map the obstacles in the robot’s workspace into its configuration space [3, 
4, 5, 13, 14]. These C-space obstacles represent those configurations of the moving object 
that would cause collisions. Free space is defined to be the complement of the C-space 
obstacles. 

Motion planning requires an explicit characterization of the robot’s free space. The 
characterization may not be complete, for example, it may cover only a subset of the 
free space. But, without a characterization of the free space, one is reduced to trial 
and error methods to find a path. In this paper we show how to compute approximate 
characterizations of the free space for simple manipulators. By simple manipulators we 
mean manipulators composed of a non-branching sequence of links connected by either 
revolute or prismatic joints (see [18] for a treatment of the kinematics of simple manip¬ 
ulators). We restrict the position of link zero of a simple manipulator to be fixed. Most 
industrial manipulators (not including parallel-jaw grippers) are simple manipulators in 
this sense. 

The C-space obstacles for a manipulator with n joints are, in general, n-dimensional 
volumes. Let C denote an n dimensional C-space obstacle for a manipulator with n joints. 
We represent approximations of C by the union of n — 1 dimensional slice projections [13, 
14]. Each n - 1 dimensional configuration in a slice projection of C represents a range 
of n dimensional configurations (differing only in the value of a single joint parameter) 
that intersects C. 

A slice projection of an n dimensional C-space obstacle is defined by a range of values 
for one of the defining parameters of the C-space and an n — 1 dimensional volume. Let 
the q = (< 7 !,..., q n ) denote a configuration, where each is a joint parameter, each of 
which measures either angular displacement (for revolute joints) or linear displacement 
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Figure 3. Slice Projection of a three-dimensional obstacle into a list of two-dimensional slices 
that are in turn represented by one-dimensional slices. 


(for prismatic joints). Let {q | qj £ [a,/?]} be the set of all configurations for which 
q } £ [a,/3] and let be a projection operator such that 

ft J (?1 ! • • • ) ?fl ) (<7l > • • • ; *7j — 1 ! Qj + 1 ! • • • 7 Qtl) 

Then, the slice projection of the obstacle C for values of q 3 € [a,0\ is 

ftj{C n {q | qj £ [a, (3}}) 

The definition of slice projection is illustrated in Figure 3. In the example above, joint j 
above is called the slice joint while the other joints are known as free joints. 

Note that a slice projection is a conservative approximation of a segment of an n 
dimensional C-space obstacle. An approximation of the full obstacle is built as the union 
of a number of n — 1 dimensional slice projections, each for a different range of values of 
the same joint parameter (Figure 3). Each of the n - 1 dimensional slice projections, in 
turn, can be approximated by the union of n - 2 dimensional slice projections and so on, 
until we have a union of one dimensional volumes, that is, linear ranges. This process is 
illustrated graphically in Figure 3. Note that the slice projection can be continued one 
more step until only zero dimensional volumes (points) remain, but this is wasteful. 

Consider a simple two-link planar manipulator whose joint parameters are qi and 
< 72 - C-space obstacles for such a manipulator are two dimensional. The one dimensional 
slice projection of a C-space obstacle C for qi € [a,/?] is some set of linear ranges {Ri} 
for <72 . The ranges must be such that if there exists a value of g 2 > call it u>, and a value 
q i € [<*,/?], call it f, for which (f,u/) £ C, then u> is in one of the Ri (Figure 3). 
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Figure 4. (a) Two link revolute manipulator and obstacles, (b) Two dimensional C-space with 
obstacles approximated by a list of one dimensional slice projections (shown dark). The initial 
and final positions of the manipulator are shown in the input space and the C-space. 


A representation of a configuration space with obstacles is illustrated in Figure 4b, 
for the two link manipulator and obstacles shown in Figure 4a. The actual configuration 
space is the surface of a torus since the top and bottom edge of the diagram coincide 
(0 = 2 n), as do the left and right edge. The obstacles are approximated as a set of qi 
ranges (shown dark) for a set of values of qi. The resolution is 2 degrees along the qi 
axis. 

If the manipulator has three links, its configuration space can be constructed as 
follows: 


1. Ignore links beyond link 1. Find the ranges of legal values of qx by considering 
rotations of link 1 around the fixed base. 

2. Sample the legal range of qi at the specified resolution. Do steps 3 through 5 for 
each of the value ranges of q\. 

3. Ignore links beyond link 2. Find the ranges of legal values of qi by considering 
rotating link 2 around the positions of joint 2 determined by the current value range 
of qx. 
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4. Sample the legal range of q 2 at the specified resolution. Do step 5 for each of these 
value ranges of q 2 ■ 

5. Find the ranges of legal values of q$ by considering rotating link 3 around the position 
of joint 3 determined by the current value ranges of qi and q 2 . 


Some sample slices from a configuration space computed in this way can be seen in Figure 
5. 

Note that the process described above is an instance of a simple recursive process: 


To compute C-space(t): 

1. Ignore links beyond link i. Find the ranges of legal values of g, by considering 
rotating link t around the positions of joint i determined by the current value ranges 
ofqi,..., qi- 1 - 

2. If * = n then stop, else sample the legal range of qi at the specified resolution. 
Compute C-space(» + 1) for each of these value ranges of qi. 


Observe that the basic computation to be done is that of determining the ranges of 
legal values for a joint parameter given ranges of values of the previous joints. This 
computation is the subject of Section 3. 

The recursive nature of the C-space computation calls for a recursive data structure 
to represent the C-space. The current implementation uses a tree whose depth is n - 1, 
where n is the number of joints, and whose branching factor is the number of intervals 
into which the legal joint parameter range for each joint is divided (Figure 6). The leaves 
of the tree are ranges of legal (or forbidden) values for the joint parameter n. Many of 
the internal nodes in the tree will have no descendants because they produce a collision 
of some link i < n. 

The main advantage of a representation method built on recursive slice projection 
is its simplicity. All operations on the representation boil down to dealing with linear 
ranges, for which very simple and efficient implementations are possible. The disadvan¬ 
tages are the loss of accuracy, and the rapid increase of storage and processing time 
with dimensionality of the C-space. Contrast this approach with one that represents the 
boundaries of the obstacles by their defining equations [4, 5]. Using the defining equations 
is cleaner and more accurate, but the algorithms for dealing with interactions between 
obstacle boundaries are very complex. I believe that the simplicity of slice projection 
outweighs its drawbacks. These drawbacks can be significantly reduced by exercising 
care in the implementation of the algorithms. 
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Figure 5. Three link revolute manipulator and obstacles. Samples of two dimensional 
slice projections used to approximate the three dimensional configuration space. The initial and 
final positions of the manipulator are shown in the input space and the C-space. 
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Figure 6. The recursive nature of the C-space leads to a recursive data structure: an n-level 
tree whose leaves represent legal ranges of configurations for the robot manipulator. 


3. Slice Projections for Polygons 

The key step in our approach is computing one dimensional slice projections of C-space 
obstacles. That is, determining the range of forbidden values of one joint parameter, 
given ranges of values for all previous joint parameters. We will illustrate how these 
ranges may be computed by considering the case of planar revolute manipulators and 
obstacles. We will first discuss this problem informally and then derive the solution from 
the equations of C-surfaces. 


3.1. A geometric view 

Assume that joint k, a revolute joint, is the free joint for a one-dimensional slice projection 
and that the previous joints are fixed at known values. Note that we assume, for now, 
that the previous joints are fixed at single values rather than ranges of values; we will 
see in Section 3.3 how to relax this restriction. We require that the configuration of the 
first k - 1 links be safe, that is, no link intersects an obstacle. This is guaranteed by the 
recursive computation we saw in Section 2. Given these assumptions, we need to find 
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Figure 7. Contact conditions for computing one dimensional slice projections: (a) Vertex of 
obstacle and edge of link (b) vertex of link and edge of obstacle. The circles indicate the path 
of the vertices as the link rotates around the specified joint. 


the ranges of values of the single joint parameter qk that are forbidden by the presence 
of objects in the workspace. 

The ranges of forbidden values for qk will be bounded by angles where link k is 
just touching an obstacle. For polygonal links moving among polygonal obstacles, the 
extremal contacts happen when a vertex of one object is in contact with an edge of 
another object. Therefore, the first step in computing the forbidden ranges for qk is to 
identify those critical values of qk for which some obstacle vertex is in contact with a link 
edge or some link vertex is in contact with an obstacle edge (Figure 7). 

The link is constrained to rotate about its joint, therefore every point on the link 
follows a circular path when the link rotates. The link vertices, in particular, are con¬ 
strained to known circular paths. The intersection of these paths with obstacle edges 
determine some of the critical values of qk, for example, B in Figure 7. As the link 
rotates, the obstacle vertices also follow known circular paths relative to the link. The 
intersection of these circles with link edges determine the remaining critical values for 
qk, for example, A in Figure 7. 

Determining whether a vertex and an edge segment can intersect requires first in¬ 
tersecting the circle traced out by the vertex and the infinite line supporting the edge 
to compute the potential intersection points. The existence of such an intersection is 
v a necessary condition for a contact between link and obstacle, but it is not sufficient. 
Three additional constraints must hold (Figure 8): in-edge constraint - the intersection 
point must be within the finite edge segment, not just the the line supporting the edge; 
orientation constraint - the orientation of the edges at the potential contact must be 
compatible, that is, the edges that define the contact vertex must both be outside of the 
contact edge; reachability constraint - for non-convex objects, there must not be other 
contacts that prevent reaching this point. 

The in-edge constraint can be tested trivially given the potential contact point and 
the endpoints of the contact edge. Since we know that the contact point is on the line 
of the edge, all that remains to be determined is whether it lies between the endpoints 
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Figure 8. Given the intersection of a vertex circle and an edge line, the following conditions 
must be met for a feasible contact: (a) The contact must be in the edge segment, contact 1 
satisfies this but l' does not (b) The edges that define the contact vertex must both be outside 
of the contact edge, contact 1 satisfies this but contact 2 does not. (c) The contact must be 
reachable, contact 1 satisfies this, but contact 3 does not (this condition is only relevant for 
non-convex objects). 


of the edge. This can be done by ensuring that the x and y coordinates of the contact 
point are within the range of x and y coordinates defined by the edge endpoints. Note 
that for contacts involving link edges and obstacle vertices, the position of the endpoints 
of the link edge must be rotated around the joint position by the computed value of the 
joint angle at the contact. 

The orientation constraint can also be tested simply. All that is required is that the 
two edges forming the contact vertex be on the outside of the contact edge. Polygon 
edges are typically oriented so that they revolve in a counterclockwise direction about 
the boundary. Therefore, the outside of the polygon is on the right of the edge as we 
traverse the boundary. Given this, the feasibility of a contact can be verified simply by 
comparing the absolute orientations of the edges involved in the contact. 

The reachability constraint, on the other hand, requires examining all the contacts 
of the link with a given obstacle that satisfy the first two constraints. For each contact 
angle q we determine whether values of q * greater than q cause collision or whether values 
less than q cause collision (Section 3.2). The contact angles together with the collision 
directions can be merged to form the ranges of forbidden values for qk■ This process is 
illustrated in Figure 9. 


3.2. Derivation using C-surfaces 

The two types of contacts (vertex-edge and edge-vertex) give rise to the two basic type 
of C-space boundary (hyper-)surfaces [3, 4, 5, 14]. One type of C-surface {type A) 
characterizes the configuration of the moving object for which a vertex of the stationary 
obstacle is in contact with the infinite line supporting an edge of the moving object. The 
other (type B ) characterizes the configuration of the moving object for which the infinite 
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Figure 9. Constructing ranges of forbidden values using the potential contact angles and the 
collision directions. 


line supporting an edge of the stationary obstacle is in contact with a vertex of the 
moving object. The equations of such surfaces are parameterized by the configuration 
parameters of the moving object. For planar polygons, x, y, 0 can be used as configuration 
parameters; for manipulators, the ft are the configuration parameters. 

For a revolute joint, choose the coordinate system to be located at the joint. The 
coordinate representation of all of the vectors will be relative to this coordinate system. 
We represent a line supporting an edge by an equation of the form: n-x+d = 0. Where n 
is the (outward pointing) unit vector that is normal to the line and d is the perpendicular 
distance to the edge from the origin. The condition for a vertex v being in contact with 
such a line is simply n • v -i- d = 0. 

For a type B contact, we are given a link vertex whose initial position vector (for 
q k = 0) is v and an obstacle edge whose line equation is n • x + d = 0. If the link angle 
is qk, the coordinates of the rotated link vertex are: 

v' = ( v x cos q k - v y sin q k , v x sin q k + v y cos q k ) 

Substituting into the plane equation yields a simple trigonometric equation in q k (all the 
other terms are constant): 

(n x v x + n y v y )cosftfe + ( n y v x - n x v y ) sing* + d — 0 (1) 

From the definition of the scalar and vector product, we have that 

n x v x + n y t; y = ||v|| cos<j>, n y v x - n x v y = ||v|| sin^ 

where <j> is the angle between n and v. From this, it is clear that the C-surface equation 
is merely 


j|v|j cos(ftt - <f>) = -d 


The solution to this equation is: 
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Figure 10. Illustration of terms used in (1) - (4) 


q k = cos 1 +<t> ( 2 ) 

Figure 10 illustrates this situation. There is one such C-surface for each combination of 
link vertex and obstacle edge. Of course, only convex vertices need be considered; no 
contact is possible at a concave vertex. 

Using the same notation, except that the edge is a link edge and the vertex an 
obstacle vertex, the equation for a type A C-surface is 

(n x v x + riyVy) cosg fc - ( n y v x - n x v y ) sinq k + d = 0 (3) 

The only difference is the sign of the coefficient of sin q k , this arises from the fact that we 
are thinking of the obstacle vertex as counter-rotating while the link stands still. That 
is, the direction of rotation of the vertex is the opposite of q k \ this changes the sign of 
the sine of the angle. The solution to this equation is: 

q k = cos -1 - 4> ( 4 ) 

There is one such C-surface for each combination of (convex) obstacle vertex and link 
edge. 

Note that there are generally two solutions to each of the equations (arising from the 
arccosine) since they correspond to intersections of a circle traced out by a vertex and an 
infinite line supporting an edge. These solutions, however, do not necessarily represent 
feasible contacts between the link and an obstacle. The remaining constraints illustrated 
in Figure 8 must also be satisfied. Of course, when the magnitude of the argument to 
the arccosine is greater than one, this indicates an infeasible contact, that is, the line is 
beyond the reach of the vertex. 

The in-edge constraint can be checked, as described before, by computing the co¬ 
ordinates of the intersection point and the positions of the edge endpoints, given the 
computed values of q k . 
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Figure 11. Testing the orientation constraint for polygonal contact 


The solutions obtained from equations (2) and (4) must also satisfy the orientation 
constraint. One way of testing this constraint is by ensuring that the polygon edges that 
intersect at the contact vertex both point outward from the contact edge. If ei and e 2 
are the edge vectors pointing away from the vertex (Figure 11), then the orientation 
constraint boils down to 


sign(n • e*) > 0, sign(n • e 2 ) > 0 
where sign(x) = x/|x| for i^O and 0 otherwise. 

The reachability constraint is handled as described in Section 3.2. To do that, we 
must be able to tell whether an increase in qk will move the contact away from the contact 
or further into contact. This can be done by computing the derivative of Eqs. (1) and 
(2). The left hand side of these equations is a measure of the perpendicular distance of 
a vertex from an edge. The sign of the derivative of this distance with respect to q * will 
indicate whether a change in q * will move further into contact or away. For example, the 
sign of the derivative of the distance for a type B contact (Eq. (1)) is determined by the 
sign of - sin(^fc + <f >) evaluated at the value of qk that gives rise to contact. 

3.3. The effect of ranges of joint angles 

Our discussion thus far has been limited to situations where all the joints except the 
last have known fixed values. The definition of one-dimensional slice projections allows 
all the joints, save one free joint, to be within a range, not just a single value. We can 
readily convert, the slice projection problem (for ranges of joint values) to the simpler 
crossection projection problem (for single joint values) we have already discussed. The 
idea is to replace the shape of the link under consideration by the area it sweeps out when 
the joints defining the slice move within their specified value ranges [13, 14]. Any safe 
placement of the expanded link represents a range of legal displacements of the original 
link within the specified joint ranges. 

In most cases, instead of computing the exact swept volumes, we can use a very 
simple approximation method. Assume the manipulator is positioned at the configuration 
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Figure 12. The manipulator link can be grown by a radius 6k' the maximum cartesian 
displacement of any point on the link in response to joint displacements e, for i < k. 


defined by the midpoint of all the joint value ranges specified for the slice projection. 
Compute the magnitude, 6k, of the largest cartesian displacement of any point on link k 
in response to any displacement within the specified range of joint values. If we “grow” 
each link by its corresponding radius 6k, the grown link includes the swept area. 

A polygonal approximation to the grown link can be obtained by computing the 
“set sum” or “Minkowski sum” of the link and a polygon enclosing a circle of radius 6 
[14]. An example of such an grown manipulator can be seen in Figure 12. 

We can illustrate this approach by considering how to compute 6k for a planar 
manipulator composed of n revolute joints. The motion of a joint affects the displacement 
of all subsequent links. Therefore, the maximum cartesian displacement of each link 
depends on the maximum total distance from any point on a link to the base joint and 
on the maximum angular displacement of the link. The maximum distance of points on 
link k is the sum of the distances between all previous joints plus the distance of any 
point on link k from joint k. The angular displacement of link A; in a planar revolute 
manipulator is also the sum of the angular displacements of all the previous joints. Given 
the distance d and the angle ff, the magnitude of the displacement (chord of a circle) is 
y/2(l - COS0). 

Let the allowed angle range for qk be ± e*; let rk be the maximum distance of 
any point on link k from joint k ; and let Ik be the distance from joint k to joint k + 1. 
The value of 6k is 
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Because the last link’s motion is never quantized when computing the C-space, we have 
that: €k = 0. This value of 8 k is very conservative; it is the largest displacement anywhere 
in the workspace. In fact, it corresponds to the displacement in a link when all the pre¬ 
vious links are fully outstretched, that is, all the ctj = 0 ,j < k. Different configurations 
would yield smaller values of Sk¬ 
in Figure 12 the relevant parameter values are: e\ — 2°, €2 = 2°, €3 = 0,/o = 0 ,/1 = 
17.0 ,12 = 17.0, r x = 18.44, f 2 = 17.26, r$ = 5.385. Therefore, the values of the 8 k are: 

= 0.644 ,82 = 2.39 ,83 — 2.749. Note the growth in the value of 8 k as the distance from 
the base increases. Because of this, one might want to choose a finer quantization for 
joints associated with long links near the base, for example, joint two in our example. 

In some applications, if the €k are small, it may be preferable to ignore the effect of 
small efc during planning and simply check the resulting path for collisions. Of course, 
if the joint ranges e* are large, these gross approximation may be too conservative and 
the exact swept volume should be used. 


3.4. Prismatic joints 

The discussion above has concentrated on revolute joints but the approach is not limited 
to them. If any of the joints are prismatic, only the computation of one dimensional 
slices will be different and, in fact, it will be simpler. 

As before, the key problem is computing the critical value of the joint parameters 
for which a link is in contact with an obstacle. These contacts involve contact of a vertex 
and an edge. So, as in the case for revolute joints, we need is determine the locus of 
motion of link vertices relative to obstacle edges and the locus of motion of obstacle 
vertices relative to link edges. For links actuated by revolute joints, we have seen that 
the vertices trace out circles. For links actuated by prismatic joints, the points on the 
link trace out lines. Potential points of contact occur where the lines defined by the 
motion of the vertices intersect the edges. This operation replaces the intersection of 
circles and lines in the preceding discussion. 

The points of intersection must still satisfy the in-edge, orientation, and reachability 
constraints. Note that the in-edge constraint must now be modified to check, not only 
that the intersection point is within the finite edge segment of the polygon, but also that 
the contact is within the range of motion of the joint. 


4. Slice Projections for Polyhedra 

The basic approach described in Section 3 carries over directly to three dimensional 
manipulators and obstacles. There is, however, one significant difference: there are three 
types of contacts possible between three dimensional polyhedra. The three contact types 
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are: (type A) vertex of obstacle and face of link, (type B) vertex of link and face of 
obstacle, and (type C) edge of link and edge of obstacle. 

Let us consider type B contacts first. Each revolute joint is characterized by an axis 
of rotation. As the joint rotates, link vertices trace circles in a plane whose normal is 
the joint axis. The intersection of this circle with the plane supporting an obstacle face 
defines two candidate points of contact (see the appendix). As in the two-dimensional 
case, possible contacts must satisfy three constraints to be feasible: in-face constraint - 
the contact must be within the obstacle face, orientation constraint - all of the link edges 
meeting at the vertex must be outside of the obstacle, and reachability constraint - for 
non-convex polyhedra, there must not be any earlier contacts that prevent reaching this 
one. 

The in-face constraint can be checked using any of the existing algorithms for testing 
whether a point is in a polygon. The orientation constraint can be enforced by checking 
that the dot products of the face normal with each of the vectors from the contact vertex 
to adjacent vertices is positive [5]. The reachability constraint is enforced exactly as in 
the two-dimensional case by merging the forbidden angle ranges. 

Type A contacts are handled analogously to type B contacts except that now the 
vertex belongs to an obstacle and the face to a link. The axis of rotation is still that of 
the manipulator joint. 

Detecting type C contacts require detecting the intersection of a line (supporting 
a link edge) rotating about the joint axis and a stationary line (supporting an obstacle 
edge). The solution for this case can be found in the appendix. Of course, an intersection 
point must be inside both edge segments to be feasible. There is also an orientation 
constraint which is a bit more difficult to derive than those for type A and B contacts 
but not particularly difficult to check (for the derivation, see [5]). 

The appendix shows the details of these computations. 


5. Free Space Representation 

Having obtained a conservative approximation of the C-space obstacles, the free space is 
simply the complement of all the obstacles. Since the obstacles are ultimately represented 
as sets of linear ranges, the complement is trivial to compute. A two dimensional free 
space, for example, will be represented as a list of one dimensional slices. Each slice 
represents the ranges of legal values of q 2 for some small range of values of q\. This is 
in itself a reasonably convenient representation of the free space but not very compact. 
If we were to try to find paths through the individual slices a great deal of time would 
be wasted searching through nearly identical slices. A more compact representation is 
called for, one that captures some of the coherence between adjacent slices. 

The free space representation used in the current implementation is made up of 
regions. A region is made up out of overlapping ranges from a set of adjacent slices 
(Figure 13). The area of common overlap of all the slices in a region is rectangular and 
called the region’s kernel. In practice, we require some minimum overlap between slices 
in the same regions to avoid very narrow kernels. 



1,3 (a) ^ (b) 

Figure 13. (a) Region definition for two link C-space. The rectangular regions are the region 
kernels. The shaded area shows region 1?2- (b) Region graph corresponding to the regions 
in part A. The link labels indicate the existence of a common boundary in the q\ and/or q 2 
directions. 

Free space regions are non-convex and so points within the region may not always be 
connectable by a straight line. There is, however, a simple method for moving between 
points within the region: move from each point along its slice to the edge of the kernel 
and connect these kernel points with a straight line. 

To search for a path between points in different regions requires representing the 
connectivity of the regions. We build a region graph where the nodes are regions and 
the links indicate regions with common boundary. Associated with each region are a 
set of links to adjacent regions, each link records the area of overlap. Regions have 
neighbors primarily in the qi direction; for these neighbors, the range of qi values at the 
common region boundary is stored with the link. By construction, regions only have q% 
neighbors at the 0 = 2x boundary, anywhere else the region is bounded above and below 
by obstacles. 

In general, each n dimensional slice is represented as a list of n — 1 dimensional 
slices and one dimensional slices are a list of ranges of joint values. We have seen that 
two dimensional regions are constructed by joining neighboring one dimensional slice- 
projections. In principle, we could construct three dimensional regions by joining neigh¬ 
boring two dimensional regions, and so on. Instead, for three dimensional C-spaces we 
simply build two dimensional regions for each range of values of the first joint parameter 
and represent the connectivity among these regions in the region graph (Figure 14). The 
connectivity is determined by detecting overlap between region kernels in neighboring 
two dimensional slices, that is, slices obtained by incrementing or decrementing the first 
joint parameter. When overlap exists, the area of overlap is associated with the corre- 
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Figure 14. Region connectivity for three dimensional slices; regions can have neighbors in qi 
direction. 


sponding link in the region graph. This method is readily extended to n dimensional 
slices by considering as neighbors slices obtained by incrementing or decrementing one 
of the first n — 2 joint parameters used to define the two dimensional slice. 

The main feature of this region representation is that it exploits the coherence of the 
free space thus, for example, it does not introduce many arbitrary divisions in the free- 
space such as are introduced by octree-type representations [7], Exploiting the natural 
coherence has a number of practical advantages. The main result is the compactness of 
the representation: very few regions are required to represent rather complex free spaces. 
Another important result is low branching in the region graph: each region has relatively 
few neighbors. These characteristics of the representation also make possible some of the 
heuristic search technique described in section 7. 

6. Searching for a Path in the Region Graph 

In this section, I describe a technique for searching a region graph, This technique applies 
to searching any subset of the C-space; it is not necessary that the complete C-space be 
examined before any searching is done. Section 7 describes some heuristic strategies for 
limiting what parts of the C-space are actually explored. 

Path searching is done by an A* search in the region graph from the region containing 
the start point to the region containing the goal point. During the search, a list of search 
nodes is kept. Each search node is associated with some intermediate region in the region 
graph and represents a set of regions connecting the start region to that intermediate 
region. For each node, we also keep track of an entry point on the region boundary that 
represents the location where the robot path would enter the region. When a search 
node is expanded by extending the region path to an adjacent region, the entry point 
is moved to the closest point on the common boundary between the two regions. The 
entry point to the next region becomes the exit point for the current region. 

To carry out the search we must associate with each search node an actual distance 
covered and an underestimate of the remaining distance to the goal. We use the distance 
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Figure 15. (a) Regions for example in Figure 4 (b) Path found between start (1) and goal (4) 
configurations (c) Some intermediate configurations. 

between entry points to define the distance between two regions and the underestimate 
is the distance between the entry point and the goal. Of course, these distances are 
based on differences between the joint parameters modulo 2x. Once having found a list 
of regions connecting the start to the goal, the actual path is obtained by connecting the 
entry points and exit points of the regions. The entry point of the start region is the 
start point and the exit point of the goal region is the goal point. 

A typical path found by the algorithm using the simple strategy described above is 
shown in Figure 15. The paths tend to be jagged; some simple postprocessing to smooth 
the path would be desirable and is currently under investigation. On the other hand, 
because of the compactness and low branching of the region representation, searching for 
a path tends to be very fast (less than half a second for two dimensional C-spaces). 

7. Heuristic Subsets of the C-space 

Having built a C-space, it may be searched repeatedly for different paths. Changes to 
the environment, however, will cause parts of the C-space to be recomputed. In rapidly 
changing environments, it may not be appropriate to compute the complete C-space 
since only small sections of the C-space will ever be traversed. This section describes 
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Figure 16. The first three links of the manipulator of Figure 1. The last three links and the end 
effector have been replaced by a simple bounding box. 


experience with a number of simple heuristic strategies that help select the subset of the 
C-space relevant to a particular path. 


7.1 Decoupling the degrees of freedom 

The path shown in Figure 1 was computed using two simple heuristics to subset the 
C-space: First plan a path for the first 3 links and a simple conservative approximation 
of the rest of the manipulator (the last three links, the end-effector and the load), see 
Figure 16. The origin and goal for this path are chosen to be the points in free space 
closest to the (projection of the) actual origin and goal. Note that these points may 
differ from the actual origin and goal in all of the joints. Having found such a path, there 
remains finding paths in the six-dimensional C-space between the actual origin (resp. 
goal) and the origin (resp. goal) of the path. For all these paths, we compute only the 
portion of the C-space bounded by the joint values of the origin and goal configurations. 

This strategy has the effect of nearly decoupling the degrees of freedom. The six¬ 
dimensional planning is confined to the areas near the origin and goal. Of course, this 
strategy will fail to find a path in the worst case, but this strategy has proven to be 
reliable and efficient in most practical situations. 


7.2 Modifying low-dimensional paths 

One alternative approach for searching for a path from configuration (q s> i,.. -,q s ,n) to 
configuration (q 9i i,..., q g , n ) is the following: 
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1. Ignore all but the first two joints. Build the free space for this reduced manipula¬ 
tor. Find a sequence of free space regions that contain a path from (?e,i,?e, 2 ) to 
(9 Sl i>?s,2)- Let * = 3. 

2. Expand the portion of the free space included in the regions found so far to incor¬ 
porate the link t. That is, quantize the range of values in the region and use them 
to compute one dimensional slice projections of the C-space for the first t links. 

3. Search for a sequence of free space regions that contain a path from (g«,i,. • • ,q s ,i) 
to ■ ■ ■, Qg,i)- If i = n then stop; else increment i and go to step 2. 


This strategy is illustrated graphically in Figure 17. The idea is to focus on rele¬ 
vant sections of the configuration space by first finding paths for successively “longer” 
manipulators, starting with a two-link manipulator and going all the way through to a 
manipulator with n links. 

It is important that at each stage we consider not just a single path of the “shorter” 
manipulator, but a sequence of regions that span all the free-space between a set of 
obstacles. The addition of link n will typically change the required path for the first 
n — 1 links. Therefore, the search space should include more than a single path so as to 
avoid the need to perform a backtracking search. 

This method has been implemented. The path for the four degree-of-freedom manip¬ 
ulator shown in Figure 18 was found by this technique. The technique leads to significant 
time savings on problems involving more than two degrees of freedom. Of course, since 
the complete configuration space is not computed, the time to plan a subsequent motion 
in the same workspace will be as long as that for the initial motion. 

8. Discussion 

The main advantages of the algorithm described here are: it is simple to implement, 
it is fast for manipulators with few degrees of freedom, it can deal with manipulators 
having many degrees of freedom including redundant manipulators, and it can deal with 
cluttered environments and non-convex polyhedral obstacles. The total wall-clock time 
to compute the C-space obstacles and then plan a path for the two-link example shown 
in Figure 4 and 15 is six seconds on a Symbolics 3600 Lisp Machine with floating-point 
operations performed in software. These times could be improved by carefully re-coding 
the algorithm, but they are already quite a bit faster than a human using an interactive 
programming system (on-line or off-line). 

The main disadvantages of the algorithm are: the approximations introduced by the 
quantization may cause the algorithm to miss legal paths in very tight environments, and 
the rapid growth in execution time with the number of robot joints. This last drawback 
is probably inherent in any general motion planner; the worst-case time bound will be 
exponential in the number of degrees of freedom [19]. 

The performance of this algorithm shows that motion planning algorithms can be fast 
enough and simple enough for practical use. I believe that in many applications automatic 
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Figure 17. It is possible to alternate between searching for a path and building the C-space. (a) 
A three degree of freedom problem. Paths for reduced manipulators (b) limit which part of the 
full-dimensional C-space (c) needs to be searched. 

motion planning will be more time effective than interactive off-line programming of 
robots. In fact, the planning times will probably be on the order of the times required 
to perform hidden surface elimination in graphics systems. 










Figure 18. (a) Initial and goal configurations for a two-dimensional manipulator with four 

degrees of freedom, (b) Path found by the algorithm in Section 7.2. 
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Figure 19. Joint coordinate system. 
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Appendix: Computing contact angles for polyhedra 

In what follows we assume that we are dealing with a convex polyhedron describing link 
k and an obstacle polyhedron (not necessarily convex). The coordinate system is chosen 
so that the origin corresponds to the position of revolute joint k and the z axis is aligned 
with the joint axis (Figure 19). The coordinate representation of all vectors is relative 
to this coordinate system. We assume that the initial position of the link polyhedron 
corresponds to q * = 0. We are interested in computing values of qk for which the link is 
in contact with the obstacle polyhedron. 

Type B contact: vertex of link and face of obstacle 

We are given a vertex of the link whose position vector is v and an obstacle face whose 
plane equation is n • x + d = 0 (n is the plane’s outward-facing unit normal). We solve 
for the angle qk that rotates the vector onto the plane. We obtain the equation for qk by 
substituting the vertex’s position, rotated by qk, into the plane equation and solving for 
<lk- 

The coordinates of the position vector for the rotated vertex are: 

v' = ( v x cos q k - Vy sin q k , v x sin q k + v y cos q k , v z ) 

Substituting into the plane equation yields n-v'+d = 0, this yields a simple trigonometric 
equation 


{n x v x + n y Vy) cos q k + ( n y v x - n x v y ) sing* = -d - n z v z (5) 

whose solution is Eq. (6). 

—l i d 1 . , , , 

qk — cos ' : + atan(n„r I - n x v y , n x v x + n y v y ) (6) 

Eq. (5) is the equation for a type B C-surface [14], Note that if we let n z = 0, then 
Eq. (5) and Eq. (1) are essentially identical. The arctangent simply computes the angle 
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between the plane normal and the projection of v on the xy plane; this magnitude is 
analogous to <f> in the planar case. Eq. (6) and Eq. (2) are also related in the same way. 

The left hand side of the Eq. (5) represents the perpendicular distance of the rotated 
vertex from the obstacle plane. The sign of the derivative of this quantity with respect 
to qk can be used to determine whether increasing or decreasing qk causes a collision. 

The orientation constraint simply requires testing whether the other endpoint of all 
the edges meeting at the contact vertex are on the outside of the plane. This is done by 
substituting the position vector of these endpoints into the left hand side of the plane 
equation and testing that the value is positive. 

Type A contact: vertex of obstacle and face of link 

We are given an obstacle vertex whose position vector is v and a link face whose plane 
equation is n • x + d = 0 (n is the plane’s outward-facing normal). The solution for qk is 
almost identical to the type A case, the only difference is the sign of the first argument 
to the arctangent. This reflects the fact that in type A contact we are treating the link 
as stationary and assuming the object is rotating in the opposite direction. This changes 
the sign of the sine of the angle. 

Type C contact: edge of obstacle and edge of link 

This case is substantially more difficult; we follow the derivation in [1]. We represent 
points on the edges parametrically in t. Therefore, points on the link’s edge are repre¬ 
sented by + v where v is the position vector of one of the endpoints of the edge and 
m is a vector along the edge (actually the difference vector between the endpoints). The 
parameter £ [0,1] parameterizes position along the edge. We can represent the vector 
along the obstacle edge similarly as t 0 n + w for t 0 £ [0,1]. 

As the edge rotates around the z axis, points on the edge trace out circles. The 
equation for points on those circles are: 

x 2 + y 2 = ( m x ti + v x ) 2 + (rriyti + v y ) 2 
z — m z t l + v z 

These can be combined by solving the second equation for ti = and substituting 

into the first to obtain: 

x 2 + y 2 = (^(z - v z ) + v x )\ (^L{z - v.) + v y ) 2 (7) 

\ m z J \ m z / 

This is an implicit equation for points on the rotation surface. 

The parametric form of the obstacle edge can be used to solve for the intersection 
of the edge with the rotation surface. 

x = n x t Q + w x , y=n y t 0 + ui y , z = n z t 0 ■+ w z 
Substituting into Eq. (7) gives a quadratic equation in t a . 
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Define the following terms: 

P = (nl + nl)ml - (m\ + mfal 

q - 2[(n x w x + n y w y )ml - (ml + my)(u>j, - - ( m x v x + m v t;j / )m 2 n z ] 

r = ( w l + w l) m l - [("»*(«>* - v z) + v x m z ) 2 + (m y (w z - v z ) + v y m z ) 2 } 

The quadratic equation that must be solved for t a is 

pt 2 0 + qt 0 + r = 0 

Having t 0 we can solve for ti since we know that the z values at contact must be equal. 
Therefore, 


n z t 0 + w z - v z 

ti = - 

m z 

Given values for and t 0 , we must first check that they are in the range [0,1] (the in-edge 
constraint), then we can compute points of intersection on each of the edges. Let 1 be 
the position vector of the intersection point on the link edge and o the position of the 
intersection point on the object edge. Then, 

qk — atan(/j,Oy l y o XJ l x o x ■f l y o y ) 

Note, however, that we have assumed, when deriving Eq. (7) that m z ^ 0. In the 
not uncommon event that it is, then all the points on the rotation surface have z = v z 
and so will the intersection point with the obstacle edge. We can use this to obtain 
t Q = ■ We can then solve for by using the fact that the contact point on the link 

edge will be on the same circle as the contact point on the obstacle edge: 

(m x t t + v x ) 2 + ( m y ti + v y ) 2 = ( n x t 0 + u;,) 2 + ( n y t a + w y ) 2 

But, we know the value of t a so this is a quadratic equation for tj. Given the values of 
ti and t 0 we can solve for qk as above. 

In addition to solving for the value of the contact angle we must compute the deriva¬ 
tive of the distance to the contact as a function of qk- As before, this will determine 
whether the contact angle is a potential upper or lower bound for a contact range. Un¬ 
fortunately, this is not as simple as it is for type A and B contacts. 

When the two edges are in contact, any motion component perpendicular to both of 
them will cause a collision while a component of motion along either edge will not cause 
a collision. The direction perpendicular to both edges is simply the cross product of the 
two edge vectors (given the link edge rotated to the contact angle): 

c(?fc) = m (q k ) x n 

The qk dependence has been indicated explicitly. 

One problem here is that we do not know whether c is outward pointing or not. We 
can decide that by dotting c with a direction known to point into the link. If ej and e 2 
are the direction vectors of edges meeting the link edge at one of its vertices (see Figure 
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Figure 20. The definition of the e, for t = 1, 2, 3, 4 


20), then ei + e 2 is a direction pointing into the link volume. Let k = sign(c • (ej + e 2 )), 
then kc is the outward-pointing normal we require (see [5] for a careful derivation). 
Given the value of k, the type C C-surface equation can be written as [14]: 

kc{q k )-(w- v(g fc )) = 0 (8) 

Differentiating the left hand side with respect to q k yields 

(( w y m x - w x m y )n z + (d z m x - v x m z )n y - ( d z m y - t> v m 2 )n I )fcsing fc + 

((tv y m y + w x m x )n z + (d z m y - v y m z )n y - ( d z m x - v x m z )n x )kcosq k 

where d z = v z - w z . If this derivative is negative then increasing q k will cause a collision. 

We are not done yet. We must guarantee that the contact satisfies the orientation 
constraint. The following are necessary and sufficient conditions for this [5]: 

s = sign(c • ei) = sign(c • e 2 ) 
s' = sign(c • e 3 ) = sign(c • e 4 ) 
s ^ s' 

These conditions are analogous to the type A and B cases. 

Although the derivation of the type C case is a bit involved, the actual amount of 
computation involved is not large. 
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